In the Claims 



Claim 1 (canceled) 

Claim 2 (previously presented) : The system estimation 
method according to claim 7, wherein the processing section 
calculates the existence condition in accordance with a 
following expression: 



i + 



- — ^HjHi > 0, i = 0,...,k (17) 



Claim 3 (previously presented) : The system estimation 
method according to claim 7, wherein the processing section 
calculates the existence condition in accordance with a 
following expression: 

-gEi + /)7/>0, i = 0, - - - , A; (18) 

.*■ * -r* 
— here, 

where the forgetting factor p and the upper limit value Yf have 
a following relation: 

0<p = l- %(Yf)^l, where %(Yf) denotes a monotonically 
damping function of Yf to satisfy %(1) = 1 and %(°°) = 0. 



Claims 4-6 (canceled) 
Claim 7 (currently amended) 




: A system estimation method, 

sound zystera or sound, field 
for making state estimation 
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robust and optimizing a forgetting factor p simultaneously in 
an estimation algorithm, in which 

for a state space model expressed by following 
expressions : 

x k+ i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 

disturbance including the system noise w k and the observation 

noise v k and is weighted with the forgetting factor p is 

suppressed to be smaller than a term corresponding to a 

previously given upper limit value Yf > and 

the system estimation method comprises: 

a step at which a processing section inputs the upper 

limit value Yf' the observation signal y k as an input of a 

filter and a value including an observation matrix H k from a 

storage section or an input section; 

a step at which the processing section determines the 

forgetting factor p 

a-eee •• ~ r e : .* •. • : • •• . . .v. : - — •••• . '. ' ; as a f o.l lowing 
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a step of executing a hyper Hoo filter at which the 
processing section reads out an initial value or a value 
including the observation matrix H k at a time from the storage 
section and obtains a filter gain K S/k by using the forgetting 
factor p and a gain matrix K k and by following expressions (20) 
to (22)~- \ 



nd (22), ; 



*k\k = + K Sik {y k - H fc x fc _i M ) (20) 

K Sik = K k (:, 1)/^(1, 1) , K k = p^KkRjJ-^J.Rl (21) 



Q(k) = 



p-\K k R~Pi l 



Rk = R\J,Rh Rl 



9* 0 



0 p5 7 /J 



Jl = 



%k\k-\ - Ek\k-lEk\k-\ 



Re,k = Rk + C k S k \k-iCl, Ck = [ ^* j , R e< k = RlkJlRjp *0|0 = £o 



■[: 



(23) 



0(k) denotes a J-unitary matrix, that is, satisfies 
0(k)J0H(k) T = J, J = (Ji © I), I denotes a unit matrix, 
Kk(:,l) denotes a column vector of a first column of 
the matrix Kv . 

-PY f \ 



here, 

x^ k | k : the estimated value of the state x k at the time k using 
the observation signals y 0 to y k , 
y k : the observation signal, 

F k : the dynamics of the system, I for simplification, 

K Sik : the filter gain, 

H k : the observation matrix, 

E A k| k : corresponding to a covariance matrix of an error of x^ k!k 
0(k): the J-unitary matrix, and 
R e , k : an auxiliary variable. 
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a step at which the processing section stores an 
estimated value of the state x k by the hyper H„ filter into the 
storage section; 

a step at which the processing section calculates an 
existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hi or 
the observation matrix Hi and the filter gain K s ,i, and 

a step at which the processing section ^ 



e •: • . . • . • ; . • : ed : t he step of 
executing the hyper Hoo filter. -r 



N-l 

dk = Y.Uk\u (k -i), /c = 0, 1,2, (34) 



is Caiiceileci by Lhe obiia i r;ed pse.jdc-ecrLO. 



Claim 8 (previously presented) : The system estimation 
method according to claim 7, wherein the step of executing the 
hyper H„o filter includes: 

a step at which the processing section calculates Z^k+iik 172 
by using the expression (22); 
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a step at which the processing section calculates the 
filter gain K s , k based on an initial condition of Z A k|k-i and an 
initial condition of C k by using the expression (21) ; 

a step at which the processing section updates a filter 
equation of the H=o filter of the expression (20); and 

a step at which the processing section repeatedly 
executes the step of calculating by using the expression (20) , 
the step of calculating by using the expression (21) and, the 
step of updating while advancing the time k. 



Claim 9 (currently amended) : A system estimation method-, 
Tor a ooo\roarrooo\uoo syst.ein or a sound oyo loar or oooad. f. ield 

) £ for making state estimation 

robust and optimizing a forgetting factor p simultaneously in 
an estimation algorithm, in which 

for a state space model expressed by following 
expressions : 
x k+ i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy 
gain which indicates a ratio of a filter error to a 
disturbance including the system noise w k and the observation 
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noise v k and is weighted with the forgetting factor p is 

suppressed to be smaller than a term corresponding to a 

previously given upper limit value y t , and 

the system estimation method comprises: 

a step at which a processing section inputs the upper 

limit value ft, the observation signal yk as an input of a 

filter and a value including an observation matrix H k from a 

storage section or an input section; 

a step at which the processing section determines the 

forgetting factor p 



;., j..c 



a step of executing a hyper H«, filter at which the 
processing section reads out an initial value or a value 
including the observation matrix H k at a time from the storage 
section and obtains a filter gain K s , k by using the forgetting 
factor p and a gain matrix K k and by following expressions: 



**|* = **-i|fc-i + K, ik {y k - H k x k . 1]k ^) 
K,, k = K k (:,l)/R e , k (l,l) , K k = ph{K k R e } : )Rl ik 



(61) 



Rl,k C k+l L k R r £ 



e(k) (63) 



here, 0(k) denotes an arbitrary J-unitary matrix, and C k 
= Cfc+iY is established, where 



R k = Rj i J l R k 2 , Rl = 



Re,k = Rk + C k S k \ k _iCfl, C k = |^ ^* j , R rk = Rl k J Y R*j., x 0 \ 0 - afc 0 



(23) 
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here, 

x^ k | k : the estimated value of the state x k at the time k using 

the observation signals y 0 to y k , 

y k : the observation signal, 

K s , k : the filter gain, 

H k : the observation matrix, 

0(k): the J-unitary matrix, and 

R e , k : an auxiliary variable. 

a step at which the processing section stores an 
estimated value of the state x k by the hyper Hoo filter into the 
storage section; 

a step at which the processing section calculates an 
existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hi or 
the observation matrix Hi and the filter gain K S;i , and 

a step at which the processing section ^ % 



the step of 

executing the hyper filter, - 



dk = Y, "hi[k\u <k -i), * = 0, 1, 2, - ■ • (34) 



Claim 10 (previously presented) : The system estimation 
method according to claim 9, wherein the step of executing the 
hyper H„ filter includes: 

a step at which the processing section calculates KT k 
based on an initial condition of R e , k +i, Rr,k+i and L~ k+ i by using 
the expression (63); 

a step at which the processing section calculates the 
filter gain K s , k based on the initial condition and by using 
the expression (62); 

a step at which the processing section updates a filter 
equation of the H„ filter of the expression (61); and 

a step at which the processing section repeatedly 
executes the step of calculating by using the expression (63), 
the step of calculating by using the expression (62), and, the 
step of updating while advancing the time k. 

Claim 11 (currently amended) : A system estimation method,. 

- v o \ - or nc.isvv compel for making state estimation 

robust and optimizing a forgetting factor p simultaneously in 
an estimation algorithm, in which 

for a state space model expressed by following 
expressions : 
x k+ i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 
w k : a system noise, 
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v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 

disturbance including the system noise w k and the observation 

noise v k and is weighted with the forgetting factor p is 

suppressed to be smaller than a term corresponding to a 

previously given upper limit value Yf# and 

the system estimation method comprises: 

a step at which a processing section inputs the upper 

limit value Yf- the observation signal y k as an input of a 

filter and a value including an observation matrix H k from a 

storage section or an input section; 

a step at which the processing section determines the 

forgetting factor p~~»<2^5*«H*K &e----*:^e---6-&&fce---si:>s-ee---«;&4e-i-----i-f v 

±^ii£tion..of__ vs.;.. 

where J( iv ; } £er^ote3 £ 
sa^isfzJLill ±J....m&.?J™l z...Sh. 

a step of executing a hyper Ho, filter at which the 
processing section reads out an initial value or a value 
including the observation matrix H k at a time from the storage 
section and obtains a filter gain K s , k by using the forgetting 
factor p and a gain matrix KT k and by following expressions: 
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x k \k = **-i|k-i + K S!k (y k - H fc xjfc_i| fc _i) (25) 
K s , k = piK k (:,l)/R e , k (l,l) (26^ 

[^^[^J-p^K-^cL (27) 
L k+1 = p- k *L k -^^R;*C kH L k (28) 

R e ,ku = ^ - C k+ iL k R r lL k C k+l (29) 
■Rr,*+i = #r,fc - LitC^R'^Ck+iLk (30) 



Where, 

C * +1 ~ [ till] ] ' = u ( k + 1 - iV )] = N* + *) u*], Hi = 0, . . . , 0] 

*.i = * + *i*ii»tf. r^J _^], i; 1|0 = diag{p 2 , , = i- x(7/) 

' 1 0 1 

~ L »= ° 0 "i* /r ° w ], tf 0 = O, = ^=p^ (31) 

here, 

y k : the observation signal, 

F :< : the dynamics of the system, * ^ c . .o • ;. cation, i 

H k : the observation matrix, 

x^ k | k : the estimated value of the state x k at the time k using 
the observation signals y 0 to y k , 

K Sik : the filter gain, obtained from the gain matrix KT k , and 
R e , k , L~ k : an auxiliary variable. 

a step at which the processing section stores an 
estimated value of the state x k by the hyper H ro filter into the 
storage section; 

a step at which the processing section calculates an 
existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hi or 
the observation matrix Hi and the filter gain K s ,i, and 

a step at which the processing section 
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leicla±...L' ' N - - v the step of 

executing the hyper H„ filter-:-.. 



dk = Y, 'hi[k\u(M), A = 0, 1, 2, - ■ • (34) 



nd an actual echo is cane si led Cy tec ocieaaeci pseudo-echo. 



Claim 12 (canceled) 



Claim 13 (previously presented) : The system estimation 
method according to claim 7, wherein an estimated value z v k!k of 
the output signal is obtained from the state estimated value 
x^ k | k at the time k by a following expression: 
z\|k = H k x^ k | k . 



Claim 14 (canceled) 



Claim 15 (currently amended) : A system estimation program 



i\U i to make state estimation 
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robust and to optimize a forgetting factor p simultaneously in 
an estimation algorithm, in which 

for a state space model expressed by following 
expressions : 

x k+ i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy 
gain which indicates a ratio of a filter error to a 
disturbance including the system noise w k and the observation 
noise v k and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value Yf > and 

the system estimation program causes the computer to 
execute : 

a step at which a processing section inputs the upper 
limit value Yf' the observation signal y k as an input of a 
filter and a value including an observation matrix H k from a 
storage section or an input section; 

a step at which the processing section determines the 
forgetting factor p •• •>* 4 -H • ^ t - - 
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a step of executing a hyper Ho» filter at which the 
processing section reads out an initial value or a value 
including the observation matrix H k at a time from the storage 
section and obtains a filter gain K S/k by using the forgetting 
factor p and a gain matrix KT k and by following expressions: 



x k \ k = afe*_ a|fc_ i + K Stk {y k - tffcXfc-iik-i) (25) 

K Stk = piK k (:,l)/R e , k (l,l) (26) 

= (27) 

L k+1 = p^L k -^^R;j t C k+l L k (28) 

R e ,k+i = Re,k ~ Ck+iL k R~lL k C k+l (29) 

i^^w^ . R^^-^[^ +1 B-^ l l ^ fc — m 

Where , 



= [ hII] J ' ** k+i = [u * +1 u( * + 1 - *)] = H k + 1) «*]> # i = Ni), o, . . . , o] 

= «i = [; _° pi} ], ii^^y n, p=i-x(7/) 

" 1 0 1 

Z„= 0 0 e fl, 0= h 1 0 1 Ro = 0, ^ = ± 0 , jfc-p**, (31) 

0 1 J L v P J 

here, 

y k : the observation signal, 

F k : the dynamics of the system, V\ ; i-^ , - uro • 
H k : the observation matrix, 

x^ k | k : the estimated value of the state x k at the time k using 
the observation signals y 0 to y k , 

K Sik : the filter gain, obtained from the gain matrix KT k , and 
R e , k , L~ k : an auxiliary variable. 
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a step at which the processing section stores an 
estimated value of the state x k by the hyper H„ filter into the 
storage section; 

a step at which the processing section calculates an 
existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hi or 
the observation matrix Hi and the filter gain K s ,i, and 

a step at which the processing section ^ 



e •: • . . • . • ; . • : ed : t he step of 
executing the hyper Hoo filter;-.:- 



d k = EM^Km, = 0, 1,2, ••- (34) 



and a.u 3Ct.'.;a.l eecec is cancelled by lice obc a 1 aed peeabo echo . 



Claim 16 (currently amended) : A computer readable 
recording medium recording a system estimation program dkc- 



c to make state estimation 
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robust and to optimize a forgetting factor p simultaneously in 
an estimation algorithm, in which 

for a state space model expressed by following 
expressions : 

x k+ i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy 
gain which indicates a ratio of a filter error to a 
disturbance including the system noise w k and the observation 
noise v k and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value Yf* and 

the computer readable recording medium recording the 
system estimation program causes the computer to execute: 

a step at which a processing section inputs the upper 
limit value Yf- the observation signal y k as an input of a 
filter and a value including an observation matrix H k from a 
storage section or an input section; 

a step at which the processing section determines the 

forgetting faCtOr p-™ ; v---: ■■■■■ ■■■■■■■■■ I...:..::..:...::....^.^^--.^--;^— ii : v: 
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a step of executing a hyper filter at which the 
processing section reads out an initial value or a value 
including the observation matrix H k at a time from the storage 
section and obtains a filter gain K S/k by using the forgetting 
factor p and a gain matrix KT k and by following expressions: 



x k \k = **-i|fc-i + K Stk (y k - H k x k -nk-i) (25) 

K Stk = p^K k (:,l)/R etk (l,l) (26) 

[^] = [^J-^fi-iLK +1 (27) 

L k+1 = P ^L k -^^R^ k C kn L k (28) 

Re,kH = R e,k ~ C k+ iL k R^J.L k ( 29 ) 

Where , 



° k+l = [ till] } ' ^ hu = u{k + l ~ N ^ = + ] ) "*]. #i = Ml)> 0. • • • - 0] 

Re^Ri + C^cl R 1 = ^ q j& 1|0 = dBg{^ , 3 ,..., „ = l-x( 7/ ) 

Z ° = J ^ W - «r fi = { "„ " N ] , K 0 = 0, * 0|0 = * 0 , = p^, (31) 
here, 

y k : the observation signal, 

F k : the dynamics of the system, 

H k : the observation matrix, 

x^ k | k : the estimated value of the state x k at the time k using 
the observation signals y 0 to y k , 

K Sik : the filter gain, obtained from the gain matrix KT k , and 
R e , k , L~ k : an auxiliary variable. 
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a step at which the processing section stores an 
estimated value of the state x k by the hyper H„ filter into the 
storage section; 

a step at which the processing section calculates an 
existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hi or 
the observation matrix Hi and the filter gain K s ,i, and 

a step at which the processing section * °-™ v 

see-? h»» ^...^ee;^..^^^^ 

the step of executing the hyper Ho» filter,™ 




dk = Y^HkWi), fc = 0,l,2,. 



(34) 



Claim 17 (currently amended) : A system estimation device, 

. v , • or • •> s conLrol, for making state estimation 

robust and optimizing a forgetting factor p simultaneously in 
an estimation algorithm, in which 

for a state space model expressed by following 
expressions : 
x k+ i = F k x k + G k w k 
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y k = H k x k + v k 
z k = H k x k 
here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

as an evaluation criterion, a maximum value of an energy 
gain which indicates a ratio of a filter error to a 
disturbance including the system noise w k and the observation 
noise v k and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value y t , and 

the system estimation device comprises: 

a processing section to execute the estimation algorithm; 

and 

a storage section to which reading and/or writing is 
performed by the processing section and which stores 
respective observed values, set values, and estimated values 
relevant to the state space model, 

further comprising: 

a means at which the processing section inputs the upper 
limit value y t , the observation signal y k as an input of a 
filter and a value including an observation matrix H k from the 
storage section or an input section; 

a means at which the processing section determines the 
forgetting factor p— * > > — v~ * i • * t • - * — .w.^..,,., % , 



19 



a means of executing a hyper H M filter at which the 
processing section reads out an initial value or a value 
including the observation matrix H k at a time from the storage 
section and obtains a filter gain K S/k by using the forgetting 
factor p and a gain matrix KT k and by following expressions: 



x k \k = **-i|fc-i + K Stk (y k - H k x k -nk-i) (25) 

K Stk = p^K k (:,l)/R etk (l,l) (26) 

[^] = [^J-^fi-iLK +1 (27) 

L k+1 = P ^L k -^^R^ k C kn L k (28) 

Re,kH = R e,k ~ C k+ iL k R^J.L k ( 29 ) 

Where , 



° k+l = [ f ] ' ^ hu = u{k + 1 " iV )] = W + 1) «*], Hi = [u(l), 0, . . . , 0] 
R e ,=R^C,S ip c], R 1 = ^ q j& 1|0 = dBg{^ , 3 ,..., , = 1- X ( 7/ ) 

Z ° = J ^ W - «r fi = { "„ " N ] , K 0 = 0, * 0|0 = * 0 , = p^, (31) 
here , 

y k : the observation signal, 

F k : the dynamics of the system, 

H k : the observation matrix, 

x^ k | k : the estimated value of the state x k at the time k using 
the observation signals y 0 to y k , 

K Sik : the filter gain, obtained from the gain matrix KT k , and 
R e , k , L~ k : an auxiliary variable. 
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a means at which the processing section stores an 
estimated value of the state x k by the hyper H„ filter into the 
storage section; 

a means at which the processing section calculates an 
existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hi or 
the observation matrix Hi and the filter gain K s ,i, and 

a means at which the processing section 



. ; • - the means of 

executing the hyper Hoo filter, -r 



if -i 

dk = Y^UQum, = 0, 1,2, ••- (34) 



and an actual echo Is cancel! eel by lice obtained pseucio- eerie . 



Claim 18 (previously presented) : The system estimation 
method according to claim 9, wherein the processing section 
calculates the existence condition in accordance with a 
following expression: 
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-0=i + pj 2 f >0, i = 0,---,k (18) 

here, 

* = 7/ ' ^ l-ffjr,/ ^ = 1 "^/) (19) 

where the forgetting factor p and the upper limit value Yf have 
a following relation: 

0<p = l- %(Yf)^l, where %(Yf) denotes a monotonically 
damping function of Yf to satisfy %(1) = 1 and %(°°) = 0. 

Claim 19 (previously presented) : The system estimation 
method according to claim 9, wherein an estimated value z v kjk of 
the output signal is obtained from the state estimated value 
x A k | k at the time k by a following expression: 
z v k | k = H k x^ k | k . 

Claim 20 (canceled) 

Claim 21 (previously presented) : The system estimation 
method according to claim 11, wherein the processing section 
calculates the existence condition in accordance with a 
following expression: 

-Q=i + ptf>0, z = 0, - - - , A: (18) 
1 2 A pHiK si , v 

where the forgetting factor p and the upper limit value Yf have 
a following relation: 

0 < p = 1 - %(Yf) ^1, where /(Yf) denotes a monotonically 
damping function of Yf to satisfy %(1) = 1 and %(°°) = 0. 



here, 

Q 
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Claim 22 (previously presented) : The system estimation 
method according to claim 11, wherein an estimated value z v kjk 
of the output signal is obtained from the state estimated 
value x^ k | k at the time k by a following expression: 
z v k | k = H k x^ k | k . 

Claim 23 (canceled) 
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